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Asymptotic properties of electromagnetic waves are studied within the context 
of Friedmann-Robertson- Walker (FRW) cosmology. Electromagnetic fields are 
considered as small perturbations on the background spacetime and Maxwell's 
equations are solved for all three cases of flat, closed and open FRW uni- 
verses. The asymptotic character of these solutions are investigated and their 
relevance to the problem of cosmological tails of electromagnetic waves is discussed. 
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1. INTRODUCTION 



It is well understood that while propagating, similar to any traveling 
wave, electromagnetic radiation interacts with the curvature of the spacetime 
[1]. Such interactions affect the propagation of these waves and cause them to 
scatter. The scattered waves, or tails, manifest themselves as partial backscat- 
tering when the source of the curvature is a localized object. However, when 
propagating throughout the universe, where the curvature exists everywhere, 
the scattering of electromagnetic waves due to the interaction with the back- 
ground curvature occurs throughout the space and at all times. The purpose 
of this paper is to understand how the background curvature of spacetime af- 
fects the propagation and properties of electromagnetic waves, particularly at 
large distances. 
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Studies of this nature have been done by Faraoni and Sonego[2,3] on 
the analysis of the propagation of scalar fields in FRW spacetime. In those 
papers, the authors have presented analytical solutions to the Klein-Gordon 
equation with non-minimal coupling and by studying the asymptotic character 
of these solutions, they have shown that standard calculations of the reflection 
and transmission coefficients are not applicable to this case. Scialom and 
Philippe [4] have also studied the asymptotic behavior of scalar fields near 
the singular points of Einstein-Klein-Gordon equation in a flat FRW universe. 
They determined the singular points of Einstein-Klein-Gordon equation and 
by analytically analyzing the asymptotic behavior of complex scalar fields, 
showed that these fields, similar to real scalar fields, support the existence of 
inflationary stage in a flat FRW universe. For vector fields, Noonan [5,6,7] 
studied the propagation of electromagnetic fields in a curved spacetime and 
showed that, while propagating in a conformally flat universe, components of 
the Faraday tensor stay tail- free. In this paper, propagation of electromagnetic 
waves in the FRW spacetime is studied and the explicit asymptotic character of 
these field in all three cases of flat, closed, and open universes will be discussed. 

To investigate the asymptotic character of electromagnetic fields one has 
to calculate the electric and magnetic components of these fields as measured 
by a standard observer, and study their properties at large distances. Here 
electromagnetic fields are considered as small perturbations on the background 
curvature, and solutions of Maxwell's equations are presented for all three 
cases of flat, closed, and open FRW universes. Asymptotic properties of elec- 
tromagnetic fields are then investigated by studying these solutions at large 
distances. 

Solutions to Maxwell's equations in a FRW universe have been presented 
by different authors [8-11]. Deng and Mannheim [12] have presented solu- 
tions to Maxwell's equations by solving these equations for spherical compo- 
nents of the electric and magnetic fields, separately. More recently, Mankin et 
al [13,14,15] have presented exact and approximate solutions to electromag- 
netic wave equations in a curved spacetime, based on the method proposed 
by Hadamard [16] and by using a higher-order Green's function for the wave 
equation. I adopt a very helpful method due to Skrotskii [17] who realized 
that electromagnetic field equations in a curved spacetime can be written in 
a non-covariant form formally equivalent to Maxwell's equations in a macro- 
scopic medium in flat spacetime. The electric and magnetic properties of 
this medium are tied to the background curvature of spacetime [17,18]. This 
method has been used by Mashhoon to solve Maxwell's equations in a closed 
expanding FRW universe [19]. 

The plan of this paper is as follows. Maxwell's equations are discussed in 
section 2. In section 3, solutions to Maxwell's equations are presented for all 
three cases of flat, closed and open FRW universes. Section 4 has to do with 
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the study of the asymptotic character of electromagnetic fields, and section 5 
concludes this study by reviewing the results and discussing their applications. 

The metric convention g a p = (+, — , — , — ) with ds 2 = g a pdx a dx^ is used 
throughout this paper. The Greek indices will indicate sums over 0,1,2 and 3 
while the Latin indices will sum over 1,2 and 3. The units used in this paper 
have been chosen such that T% = c = 1 where c is the speed of light. 



2. MAXWELL'S EQUATIONS IN FRW SPACETIME 

In a curved spacetime with metric g a p, the source- free Maxwell's equa- 
tions are given by 

[(-g) l/2 F^]^ = 0, (1) 

and 

Fa/3, 7 + F jay p + Fp Jy0l = , (2) 

where F a p is the electromagnetic field tensor, g = det(g al 3), and (,) denotes 
an ordinary differentiation. In accordance with Skrotskii's formalism [17], a 
background inertial frame is introduced with Cartesian coordinates in which 
the electric and magnetic fields are defined using the decompositions F a p — > 
(E,B), and ^gF^ -> (—3, H). That is, 

E a = F a0 , D a = (-g) 1 / 2 F 0a , (3) 

and 

B a = l -e ahc F hc , H a = ^e abc (-g)*F bc , (4) 

where e a \, c is the three-dimensional Levi-Civita symbol. Using these quantities, 
equations (1) and (2) can be written as [18] 

D a , a = B a , a = , (5) 

and 

— F> a ,0 + ^abcH c ,b = B aj0 + € a bcE c ,b = . (6) 

Equations (5) and (6) are formally equivalent to the electromagnetic field 
equations in a medium in flat spacetime with a dielectric constant e a b and a 
permutivity fx a b given by 

ab 

Sab = fiab = -(-9) l/2 — - (7) 

9 00 
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In this auxiliary medium and with the assumed Cartesian coordinate system, 
(E, D) represent the electric fields and (£> , H) are their corresponding magnetic 
fields. These vector fields are related via constitutive relations [19] 

D a = e ab E b - (G x H) a , (8) 
B a = VabH b + (G x E) a , (9) 

where G a = -g 0a /g 00 . 

The above-mentioned formalism is applicable to any curved spacetime. In 
this study the curved spacetime of interest is FRW with a metric given by 



(ds) 2 = {dtf - s 2 {t)Rl | (l ( ^ 2) + n\def + n 2 sin 2 ^(#) 2 | . (10) 

In equation (10), Rq is the radius of the model universe at some epoch to. 
The expansion parameter of this model universe is given by S(t) such that 
S(to) = 1. The product S(t)Ro in equation (10) is the cosmic scale factor, and 
k = —1,0,-1-1 corresponding to open, flat, and closed universes, respectively. 
For k = +1, 1Z ranges from to 1. 

Introducing x\ ,X2 and X3 given by r sin# cos </>, r sin 6 sin and rcos#, 
respectively, the line element (10) can be written as 

(ds) 2 =C 2 ( V ){(d V ) 2 - [f k (r)] 2 (S ww ,dx w dx w ')Y (11) 

In this equation, 77, the conformal time, is given by dt = S(t)drj, and C(rj) = 
S(t). The function fk(r) in equation (11) is equal to 



fk(r) = 



1 + k 



2R 



(12) 



where r is defined via 



Ro 



(13) 



In terms of the new coordinates (77, x\, X2, x%), the dielectric constant e ab 
and permutivity [i ab are given by 



£ab = Vab = fk(r)Sab ■ 



(14) 
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Also in this new coordinate system, G a = 0. Substituting these quantities in 
equations (5) and (6), one can show that Maxwell's equations can be written 
as [12] 



- dl 
-iVxl=f k (r)—, 



and 



fk{r)X 



= 0, 



(15) 
(16) 



where T = E + iH. Equations (15) and (16) represent the electromagnetic 
field equations for the line element (11). It is evident from these equations 
that for the metric (11), electromagnetic fields are independent of S(t). This 
is due to the explicit conformal invariance of the formalism above, and im- 
plies that in a universe with a background curvature given by equation (11), 
the physical parameters of electromagnetic fields vary adiabatically with the 
changing radius of the universe. Among the solutions of equations (15) and 
(16), only those that represent E and H at a conformal time r/, corresponding 
to t = f, are solutions to the Maxwell's equations in a FRW universe at time 
f. 



3. SOLUTIONS OF MAXWELL'S EQUATIONS 

The isotropy and homogeneity of FRW spacetime implies that it is possi- 
ble to find solutions to equations (15) and (16) for a definite angular momen- 
tum J, and its component along the X3-axis, M. Denoting such solutions by 
ZjMa k (r, rj), one can write 



-io k r\ 



(17) 



A=e,m,0 



where IjMa k ■> V) have been expanded in terms of the vector spherical harmon- 
ics Yj^i(9, 0) [19,20]. In equation (17), a k is a positive constant and A = e, m, 
representing different types of vector spherical harmonics [20]. For a defi- 
nite parity and angular momentum, the electric and magnetic components of 
XjMa k {r, i]) can be written as [21] 



and 



EjZa k (r) = Hf Mak (r ) = l£ k (r)Y^(e, 0) , 



(18) 
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where from equations (15) and (16), ij^ (r) are given by 



(Tkfk{r) \dr r J 



(20) 



(Tkrfk(r) 



1/2 



j(j+i) xr fe w, 



(21) 



d 2 2 d 
dr 2 r dr 



1 df k (r) fl , d 



/^(r) dr \r dr 



J(J + 1) 



^(0 = 



(22) 



The most general solution to equations (15) and (16) is a superposition of the 
fields (18) and (19). That is, 



J,M,crfc A'=e,m 



%l t (^?)^M CTfc (0 + (r/)4^ fe (r) 



A'* L 1 ! ?A'* 



(23) 



and 



^) = E 



E 



%M CTfe (vW Mak (r) + q^; fe (r?)4^ (r) 



rA' 



J,M,<Jk \'—e.m L 



(24) 



where * indicates complex conjugation and from equation (15), aj M(J {rj) satisfies 



(25) 



In the following, solutions of the differential equations (20) to (22) are 
presented for all three cases of flat, closed and open FRW universes. 
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3-a. Flat Universe 

In a flat universe where k = 0, equation (22) is written as 



cT 2 / d\ 
dr 2 r \drJ 



J(J+1) , 2 
12 + 



2™ (r) = 



(26) 



The general solution of equation (26) can be written in terms of Hankel func- 
tions h^' 2 \ That is, 



2£(r) = Ah$\ra ) + Bh^{ra ), 



(2), 



JCTQ 



(27) 



where and B are constant quantities. In equation (27), 



h ( j\ra ) = jj(ra ) + inj(ra ) 



(28) 



and 



hy(ra ) =jj(ra ) -inj(ra ), 



(29) 



where jj and nj are spherical Bessel and Neumann functions, respectively, and 
J = 0, 1, 2, 3, ... . Because nj(rcro) is singular at r = 0, solutions of equation 
(26) are proportional only to spherical Bessel functions. Using 



roc 



y vw , 



(30) 



these solutions can be written as 



X T (r) = 



2a (2J + l) 



7T 



1/2 



jj{ra Q ) , 



(31) 



where, 2™ (rao) has been normalized such that 



f 

Jo 



(32) 
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Figure 1. From top to bottom, graphs of (r) for A = 171, e, 0, and for different values 

of J and (To- The solid line corresponds to Co = 1, the dashed line corresponds to (To = 2 
and the dotted line corresponds to <Tn = 3. Because 1° (r) = for J = 0, the graphs of 

JCTQ >■ ' 

this function have been plotted for J > 1. It is also necessary to mention that the graphs of 
(A = m, J = 0) , (A = e, J = 1) and (A = 0, J = 1) meet their corresponding vertical 
axes at values equal to 5, -8 and -9, respectively. Note the different scales on the vertical axes. 



Substituting for 1™ (r) from equation (31) in equations (20) and (21), one can 
show that 



2(T 


1/2 r 


_vr(2J + 1)_ 





Jj J+1 (ra ) - (J+ tfjj.^rao) 



(33) 



and 



2£» = 



2J(J + 1)(2J + 1) 
a 7r 



1/2 



3 j ( r(T o) 



(34) 



Equations (31), (33) and (34) represent the r-dependence of the electric and 
magnetic fields given by equations (18) and (19). Figure 1 shows these func- 
tions for different values of J and ctq- From equation (34) one can see that 



7C/4 7l/2 3JT/4 % 

1 



Figure 2. Graph of J(J + l)/sin 2 X against X for J = 1, 2, 3, 4. 

2"2 (r) vanishes when J = 0. That means, Eq M(7o (f) and H^ /IcrQ (r ) will be 
pure e-type for all values of do- One can also see from figure 1 that for any 
state of A, when J is constant, as expected, increasing <7o results in increasing 
the frequency of the solutions of equation (26). This figure also shows that 
among different states of A, Ij^ (r) approaches zero more rapidly implying that 

X m (r) and X e (r) will have dominating effects on the behavior of E and H 
fields at large distances. 

3-b. Closed Universe 

In a closed FRW universe, k = 1, and therefore fk(f) = fi( r ) = 
1/{1 + [r/(2i? )] 2 }. Equation (22), in this case, is given by 

(X) = , (35) 

where uj\ = Roa±, \ = 2tan~ 1 [r/(2i? )], and iff™ (x) = i?o rl™ (r). In terms 
of the variable x, equation (35) resembles a Schrodinger equation with a po- 
tential function given by J (J + l)/sin 2 x, and eigenvalues equal to Figure 
2 shows the graph of this potential function for J = 1,2,3,4. It is evident 
from this figure that for all values of J, the solutions of equation (35) must be 
finite for all values of x £ [0, 7r]. Using this boundary condition, the solution 
of equation(35) can be written as 



d 2 

dx 2 J " 1 



(X) + 



j(j + i; 

sin 2 x 
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KM) = 2 J+1 

where J = 1,2,3,... and u>i = 2,3,4,... representing frequencies of different 
modes of electromagnetic fields for the line element (11) [8,9,19]. In equation 
(36), Cy(Q) are Gegenbauer polynomials [22] which are orthogonal over the 
interval (-1, 1) with the weight function W(G) = (1 - 6 2 )«-i/2 j v > -\/2. 

Orthogonality of (cos x) requires VP™ (x) to be normalized such that 

/ V (x)*™ (X)dx = 2™^^ . (37) 
J o 1 



(wi - j - r 

K + J)! 



. ■ L /- a 



(sin X ) J+1 C'^ (cos x), (36) 



Using 



dG 



C£(6) = 2uC^I i 1 (0), 



(38) 



and the recurrence relation 



2u(l - G 2 ) C^t^O) = (T + 2v - 1)C^_ 1 (0) - TOC^(e) , (39) 



and also from the definition of VP™ (%), and equations (20) and (21), the r- 
dependence of the electric and magnetic fields in a closed FRW universe are 
given by 



(wi + J)! 



1/2 



(1 + cos x) sin J x (cos x) , (40) 



T L ix) = -2- 



J!_ 

# 



(cji - J - 1)! 
(wi + J)! 



1/2 



(1 + cosx) sin J x 



ui cosxC wi J i (cosx) + (wi + «/)C^_ 7 _ 2 (cosx) 



(41) 



(*) = -2' 



J!_ 

i? 



J(J + l)(wi- J- 1)! 



1/2 



(wi + J)! 

(1 + cosx) sin J x C 



— J- 



^cos x)- 



(42) 
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Figure 3. From top to bottom, graphs of Rq (r) for A = 171, e, and for different values 
of J and U)\. When A = m, 0, the solid line corresponds to U)\ = .7+1, the dashed line 
correspond to Wi = J + 2 and the dotted line corresponds to Wi = J + 3. For A = e, the 
value of a>i is equal to J + 2, ,7 + 3 and J" + 4 for the solid line, dashed line and the dotted 
line, respectively. Note the different scales on the vertical axes. 



Figure 3 shows the quantity R\X^^ (r) for different values of J and uj\ . It 

is important to mention that because the index (u\ — J ± 1) in C'l +1 _ (cos x) 
can only be equal to zero or positive integers, the graphs of figure 3 contain 
degeneracies on J and u>±. As shown in this figure, for a given value of i?o, 
X (r) vanishes at \ = and 7r. This is quite expected since from figure 
2, the potential function of the Schrodinger-type equation (35) approaches 
infinity at these two values of x- This indicates that in a model closed FRW 
universe, different modes of electromagnetic radiations, which correspond to 
different values of u>i, are extended over the entire volume of the universe. 
Each of these modes is a member of a complete set of eigenfunctions with 
amplitude that increase by increasing the frequency. An electromagnetic field 
is a superposition of these eigenfunctions. 
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Figure 4. Graph of J(J + 1)/ sinh 2 It against it for J = 1, 2, 3, 4. 

3-c . Open Universe 

In an open FRW universe where k = — 1, and fk(r) = f-i(r) = 
1/{1 — [r / (2R )] 2 } , equation (22) can be written as 



du 2 



C _ («) + 



J(J + 1) 
sinh 2 it 



$ m (it) = 



(43) 



In this equation, u>_ 1 = R$(J_ 1 , u = 2tanh 1 [r/(2i? )] and <E>™ (u) = 
R rl™ (r). Equation (43) is a Schrodinger-type differential equation with 

a potential function of J(J + l)/sinh 2 -u. Figure 4 shows this function for dif- 
ferent values of J. 

With appropriate change of variables, it is possible to show that equation 
(43) is, in fact, the differential equation of a hypergeometric function. Let 



(44) 



and 



£ = - (1 — cothit) 



(45) 



Substituting for <E>™ (u) and u from equations (44) and (45) in equation (43), 
one can show that this equation can be written as 
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1 (0 

V A ^ rf£ (46) 

J(J+1) + «,,_, (1-iwjlC (6=0. 



Introducing p = J + l — iu>_ 1 , q = —J — iu>_ 1 , and s = l — iu_ 1 equation (46) 
represents the differential equation of a hypergeometric function [22] with the 
general form of 



£(i-0 g + [*-(p + q + i)t]— ^— -p?C_ 1 (0 = o. (47) 

Equation (47), as the general differential equation of a hypergeometric 
function, has three regular singular points at £ = 0, 1 and oo. In the neighbor- 
hood of each of these singularities, equation (47) has two independent solutions 
[22]. Since for all values of < u < oo, from equation (45), — oo < £ < 0, 
among all these solutions, only the ones in the neighborhood of £ = are 
valid solutions of equation (47). Within this range, equation (47) has only one 
regular singular point at £ = 0, and its general solution is given by [22] 

C , (£) = <?, s, + T>e-F(p -8 + l,q-8 + l,2- 8 ,0. (48) 

In this equation, 

^, 9 , s ,, ) = 1 + ? i + ? <E+^f + ..., (49) 

is the hypergeometric function, defined within its circle of convergence |£| < 1, 
and B and V are constant quantities. It is necessary to mention that solution 
(48) converges conditionally at £ = — 1, since Rel(s — p — q) = [23]. 

To determine the values of the two constant quantities B and T>, it is more 
convenient to study the asymptotic character of <E>™ (u) at u — > and oo. 

Using the identity 

F(p, q, s, = (1 - S ~ P ~ q Hs -p,s-q,s,0, (50) 

and considering that, as shown by equation (49), T{p, q, s, £) is symmetric on 
p and q, equation (48) can be written as 



Q (0 =B(1- £)-- F(-J, J +1,1- o 

(51) 



13 



From equations (44) and (45), one can show that equation (51) simplifies to 

—5-) " ^(-J.J + M-iw,,*) 

2 7 (52) 

+ " ^(-J, J +1,1 + 10^,0- 

As shown in figure 4, the potential function of the Schrodinger-type equation 
(43) approaches infinity when u — > 0, and it tends to zero when u — > oo. That 
implies, the solutions of equation (43) are finite for all values of it G (0, oo). The 
general solution should therefore be a superposition of ingoing and outgoing 
radiation as in equation (52). However, the physical interpretation of the 
FRW universes based on the big bang model implies that all the radiation 
must be outgoing. That means, when u — > oo, the only possible solution 
for equation (43) is an out-going wave [2,24]. This requires no in-coming 
radiation in equation (52). However, when u — > oo, the quantity £ tends to 
zero, T{p, q,s,£) — > 1, and the second term of equation (52) represents a pure 
in-coming wave. The condition above implies that, V = 0, and as a result, the 
solution of equation (43) is written as 

C = ) " ^(- J > J + 1,1-^,0- (53) 

Setting B = 1, equation (53) represents an out-going radiation with frequency 
lu_ 1 . From the definition of <E>™ (u), and from equations (20) and (21), the 

r-dependent terms of electromagnetic fields are now given by 



XI i (r) = - JL X{r) | ?(-J, J + l,l-ta;_ 1J £(r)) 

1 (55) 

+ fo J ( J + 1 ) [/_ 1 (r)]"V(-J + 1,7 + 2,2-^,^)) , 
2r^ w_ 1 (1 — tuj_ 1 ) V / 



and 



(56) 

^^(-J.J+M-zuL^r)). 
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Figure 5. From top to bottom, graphs of Rq\X^ (r)| for A = TYl, e, 0, and for different 

values of J and 0J—\. The solid line corresponds to U!—i = 1, the dashed line corresponds to 
U)_ 1 = 2 and the dotted line corresponds to U)_ 1 = 3. Note different scales on vertical axes. 



In equations (54), (55), and (56), 



X(r) = 



1 + s(r) 

1 - q(r) 



(57) 



q(r) = (v/2Rq), and the derivatives of the hypergeometric function 
T(p, q, s,£(r)) have been replaced by 



— T(p, q, s,Z(r)) = ™r(p+l, q + 1, s + 1 , ^(r)) 

Figure 5 shows the product of the absolute values of complex functions (r) 

and Rq in terms of q, and for different values of J and It is necessary to 
mention that in an open FRW universe, 1Z G (0, oo), which from equation (13) 
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dm 

dr 



(58) 




Figure 6. Graphs of Rq\X^ (r)| for A = m (left column), A = e (right column), and for 

different values of J and CJ_i. The solid line corresponds to J = 0, the dashed line corresponds 
to J = 1 and the dotted line corresponds to J = 5. As mentioned in the text, the quantities 
Rq \L™ (2-Ro)| an d -^ol-^jL (^-^o)| are independent of the value of J and are equal to 

0.5 exp(o;_ 1 7r/2). For UJ_ 1 = 1 this value is ~ 2.4 (top graphs), for 0J_ 1 = 2, it is ~ 11.57 
(middle graphs), and when U)_ 1 = 3, it attains the value of ~ 55.66 (bottom graphs). Note 
different scales on vertical axes. 
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implies < r < 2R . As shown here, for a given value of the radius of 
universe, (r) approaches zero for A = 0, and a non-zero constant value for 

A = e,m. Using identities X Y = exp(FLnX), and ±i = exp[±z (tt/2 + 2nl)], 
and choosing the principal branch on the complex plane, one can show that 
this constant value is equal to 0.5 exp(o;_ 1 7r/2) as demonstrated in figure 6. 

4. Asymptotic Character of E(r, rj) and H(r, rj) 

To study the asymptotic behavior of electromagnetic fields, one has to 
evaluate the magnitude of E and H for large values of 71. To do so, imagine 
a set of fundamental observers in FRW universes, at rest in space with con- 
stant spatial coordinates x±, x 2 and x% as given by equation (11) . At these 
coordinates, the measured components of the Faraday tensor are given by 

F (pX „) = F aP A a {p) A\ u) , (59) 

where 

A>) = (c-\v),0,0,o), (60) 

A Q (i) = (o.C-^iy^CrJ.O.o), (61) 

A Q ( 2 ) = (O.O.C-H^W.O), (62) 

A Q (3) = (O.O.O.C-WW). (63) 

are the orthonormal tetrad frames associated with these observers calculated 
from equation 

g ap = A« {p)A P {u) . ( 64 ) 



In equation (64), i]^^ is the metric of the Minkowski spacetime and g a/3 
represents the line element of FRW universes as given by equation (11). Using 
equations (3) and (4) and also equations (7) to (9), the components of Faraday 
tensor F^u) can be written as 



F {pXu) = C- 2 (rj) [f k (r)] 



-i 



/ 

-E x 

—E2 
\-E 3 



Et 


-H 3 
H 2 



E 2 
H 3 




E 3 \ 
/ 



(65) 



17 



In equation (65), (Ei,E 2 , E 3 ) and (Hi,H2, H s ) represent the components 
of the electric and magnetic fields E(f,rj) and H(r, 77), as given by equations 
(23) and (24), along x\- x 2 - and X3-axis. To study the asymptotic behavior of 
E and H, one has to study the properties of these components for large values 
of 71. Such a study is meaningful only within the context of a flat and an open 
FRW universe. In a closed FRW universe, as shown in section 3-b, at any 
given epoch, the radial parts of the solution of equations (15) and (16) extend 
within the entire volume of the universe. This characteristic of X x (r) makes 

the study of the asymptotic behavior of electromagnetic fields in a closed FRW 
universe implausible. 

From equations (18) and (19), and for given values of J,M,at and A, 
the r-dependence of electromagnetic fields is given by Xj a (r). As a result, 
the measured Faraday tensor F/ p \r v \ will have an r-dependence of the form 
Xj ak (r)/ fkix)- In a flat FRW universe, TZ = t/Rq , fo(r) = 1, and from equa- 
tions (31), (33) and (34), for a constant value of cr , the functions Tj ak (r) are 
proportional to spherical Bessel functions jj(r). For large values of 1Z, the 
asymptotic values of these functions are given by [22,25] 

jjir') ~ ^-sin(r'- (66) 

Equation (66) indicates that, as expected, in a flat FRW universe, electro- 
magnetic fields approach zero at large distances with amplitudes decreasing 
as 1Z~ 1 . A result that is also seen from figure 1. In an open FRW, as indi- 
cated by equations (54), (55) and (56), (r) approach constant values as 

r — > 2Rq. Therefore from equation (65), the measured fields approach zero 
since [/_i(r)] = {1 — [ r /(2-Ro)] 2 } —> 0. It follows that, similar to a flat uni- 
verse, in an open FRW universe also, measured fields approach zero as 
as 1Z — > 00. 

5. SUMMARY 

The results of a study of the asymptotic properties of electromagnetic 
waves in a FRW universe were presented. Electromagnetic fields were consid- 
ered as small perturbations on the background curvature of FRW spacetime, 
and the metric of the universe was written in form of a line element that, at a 
conformal time corresponding to a given epoch, is electromagnetically equiv- 
alent to FRW metric. The solutions of Maxwell's equations were obtained for 
all three cases of flat, closed and open FRW universes. 

The isotropy and homogeneity of FRW metric allows writing the solutions 
of the electromagnetic field equations in terms of vector spherical harmonics. 
Expansions to generalized spherical harmonics have also been presented by 
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Mankin et al [13-15,26], and by Laas et al [27] to construct exact and approxi- 
mate solutions to electromagnetic wave equations in a curved spacetime using 
higher-order Green's function method. 

It was shown that using appropriate transformations, the equations gov- 
erning the r-dependence of the electromagnetic fields for the type-m vector 
spherical harmonics could be written in form of a one-dimensional Schrodinger- 
type equation whose eigenvalues represent the frequencies of different modes 
of electromagnetic radiations. In a closed FRW universe, the solutions of the 
above-mentioned equation resemble standing waves, extended over the entire 
volume of the universe, implying that the study of their asymptotic character 
would be implausible. 

The asymptotic values of electromagnetic fields in a flat and an open 
FRW universe were obtained by calculating these fields in the local frame 
of an observer at large distances. Analysis of the results indicated that, as 
expected, these fields tend to zero as 1Z~ 1 at large distances. 

An application of the results of the analysis presented here would be in 
the study of the tails of electromagnetic waves in FRW universes. As shown 
by Noonan [5,6,7], the vector potential of an electromagnetic wave has a non- 
zero tail when propagating in a curved spacetime. However, in conformally flat 
universes such as FRW, because of the conformal invariance of Maxwell's equa- 
tions, the components of Faraday's tensor have no tails [5,6]. Such a conclusion 
can trivially be made in a flat universe where the background spatial curvature 
is non-existent, and can be attributed to the fact that in a (3+1) flat universe, 
the four dimensional Green's function of the spacetime has a delta-function 
character, and as a result, only motions along the light cone are supported. 
For a closed FRW universe, as mentioned in section 3-b, because the solutions 
of its corresponding Schrodinger-type differential equation resemble standing 
waves, the concept of tails is not applicable. However, such studies are quite 
plausible in an open FRW universe. For a specific electromagnetic wave, in 
this case, the tail-free nature of the components of Faraday tensor can be ex- 
amined by expanding these components in terms of EjM ak if ) and HjM ak if) as 
given by equations (18) and (19), and calculating their asymptotic values for 
large values of 71. For given values of J , M and <7o, the r-dependence of these 
fields are given by Xj^ (r)//_i(r). As shown in section 4, for large values of 

TZ, these quantities approach zero as TZ^ 1 , indicating the tail- free nature of 
electromagnetic waves. 
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